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ABSTRACT 

We calculate the scattering effects associated with the interaction of a surface gravity or /mode 
with a thin magnetic flux tube embedded in a realistically stratified medium. We find that the 
dominant scattered wave is an / mode with amplitude and phase of 1.17% and around 49° relative 
to the incident wave, compared to the values of 0.13% and 40° estimated from observations. The 
extent of scattering into high-order acoustic p modes is too weak to be accurately characterized. 
We recover the result that the degree of scattering is enhanced as (a) the frequency of the incident 
wave increases and (b) the flux tube becomes magnetically dominated. 

Subject headings: Sun: helioseismology — Sun: interior — Sun: oscillations — waves — hydrodynamics 



Introduction 



A problem of considerable interest and more recently controversy relates to the influence of magnetic 
fields on acoustic waves in the near-surface regions of the Sun. In t hese sub-photospheric magnetic regions, the 
ratio of magnetic to gas pressure could be very close to unity (e.g. Mcintosh et alj 20031 Mcintosh fc Jefferiesl 



(e.g. 



2006 ), leading to the contention that magnetic fie ld effects are systematic and significant. Mode conversion 



Spruit fc Bogdan|[l99ilBarnes fc CallvfeoOOri . is a phenomenon commonly invoked to describe magneto- 



hydrodynamic (MHD) wave interactions; from one acoustic mode to another, from acoustic modes to Alfven 
waves and so on, it is estimate d that acous tic energy is redistributed, contributing perhaps to p-modc 
absorption observed in sunspots (Braun 19951 ). 



Thin flux tubes are importa nt building bloc ks of solar magnetic activit y, intimately tied to studies of the 
solar photospheric dynamo (e.g. CattaneJ 1999 ) and flux emergence (e.g. Cheung et al.l 2006) . One way of 



deconstructing the structure and dynamics of thin flux tubes is through experimental analyses of changes in 
the statistics of t he wave field in the vicinity of the s e mag netic features. Applying techniques of time-distance 
helioseismology ( Duvall et al. 19931 ). Duvall et al. ( 20061 ) estimated the phase and amplitude of the scattered 
/ modes from observations of thousands of isolated, small magnetic features. These parameters contain the 
scattering properties of an 'average' flux tube, giving us insights into the nature of these flux tubes. 
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However, one must apply caution when invoking causative mechanisms to explain thes e statistics, taking 
care t o understand the nature of wave scatterin g in the presence of m agnetic fields (e.g. ISpruit fc Bogdanl 
19921 iBoedan fc Callvlfliil iBogdan et al.lll996t Barnes fc Callvll2000h . Once a sound theoretical basis is 
formed, observations characterizing the extent of scatter can be used to place constraints on theories of the 
structure and dynamics of flux tubes. 

In an attemp t to construct an an alytical framework to decipher the nature of wave interaction with 
magnetic re gions, IBoedan et alj (|1996T ) model a magnetic flux concentration in the thin tube limit (e.g. 
Spruitlll99lh embedded in a truncated polytrope. The flux tube is assumed to be thin enough that variations 
across its interior are neglected; moreover, the tube must be a sub- wavelength feature and the radius small 
in comparison to the pressure scale height. The parameters of the polytrope are selected so as to mimic 
the solar sub-photospheric layers as closely as possible and is truncated a little below the photosphere 
because the rapidly diminishing density makes the thin tube approximation inapplicable very close to the 
surface. These simplifying assumptions reduce the generality of the results presented here and in fact, direct 
numerical solvers might appear to be a more attractive proposition. However, while numerical simulations 
can undoubtedly be employed to address a broader range of questions, it is also important to develop an 
analytical treatment of a comparable situation, as we attempt to do in this paper. Such results serve not 
only as benchmarks for direct computations of these phenomena but also assist us in developing a theoretical 
appreciation of MHD-wave interactions. 



Bogdan et al. ( 19961 ) have demonstrated that the dominant flux tube mode is the kink mod e and that 
the / mode couples strongly with the flux tube. Moreover, the measurements of iDuvall et alJ (|2006l ) are 
for / modes. Consequently, the focus of this paper will be to quantify the interaction between an incoming 
/ mode and the resulting kink mode exhibited by the flux tube. In § [3J we describe the model in greater 
detail, followed by a discussion of the method we apply to extract the scattering coefficients in § [3] The 
theoretical values o f the phase and amplitude of the scattered / mode are compared with those obtained by 
Duvall et al. ( 20061 ) (also see appendix[B|). We also demonstrate that the scattering process is predominantly 
restricted to /-/ while the f-p conversions are very weak. We conclude with a summary of our calculations 
and a discussion of the relevance and importance of these results in § HI 



The model 



The background structure in this calculation, adapted from IBoedan et al.l (|1996T ). is an adiabatically 



stratified, truncated polytrope with index m — 1.5, gravity g = —2.775 x 10 4 cm s~ 2 z, reference pressure 
Po = 1.21 x 10 5 g cm -1 s -2 , and reference density po = 2.78 x 10 -7 g cm~ 3 , such that the pressure and 
density variations are given by, 



p( z ) =Po\ 



~0 



m+1 



and 



p( z ) ~ Po ( 



(1) 



(2) 



We utilize a right-handed cylindrical co-ordinate system in our calculations, with coordinates x = (r, 9, z) 
and corresponding unit vectors (r, 8, z). The photospheric leve l of the background mo del is at z = 0, with the 
upper boundary placed at a depth of z — 392 km. Following Barnes &; Callv ( 2000[ ), we introduce a lower 
boundary at a depth of 98 Mm. The displacement potential ^(x, t) describing the oscillation modes [t is time) 
is required to enforce zero Lagrangian pressure perturbation boundary conditions at both boundaries. This 
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upper boundary co ndition is reflective in n ature and therefore, possibly not very realistic. In comparison, 
the formalism that ICrouch fc Callvl (|l999h adopt in treating upward propagating waves at the boundary 
appears to be a more fitting description but is rather difficult to implement. While the choice of 'good' 
boundary conditions is an important issue, that of incorporating the effects of the ubi quitous magnetic field 
that pervades t he atmospheric layers is an equally, if not more important aspect (e.g. Crouch fc Callv 1999t 
Hanasogd 120071 ) . Moreover, the thin flux tube approximation rapidly breaks down as one proceeds higher 
into the atmosphere. As a first step, we have therefore chosen to (1) truncate the background model a little 
below the photosphere and (2) adopt simplistic boundary conditions, thereby not having to deal with the 
uncomfortably complex phy sics of the atmosphere. The inco ming /-mode, a plane wave, which expanded in 
cylindrical coordinates (e.g. Bogdan 1989 : Gizon et al. 2006h has a displacement eigenf unction, 'fine, of the 
form: 



(3) 



where, 



%{Kl;s) = s- l ' 2 -> i N n 



In equations §3§ and ([3]), we have introduced a slew of new symbols, 



(4) 



m — 1 



mus 2 Zo 



2k P 1 Z ' 



(5) 



u> the angular frequency of oscillation, s = —z/zg, J m (w), the Besse l function of order m and argument w and 
M K ^(w), the Whittaker function (e.g. Whittaker fc Watsonl 1980 ) with indices k, // and argument w. The 
eigenvalue K v n > and constant (^ characterizing the mode are obtained through the procedure described 
in appendix [Al The n = mode corresponds to the surface gravity or / mode, while n > represents the 
acoustic p n mode. The term N n is the normalization constant for the mode, defined as 



CP M v + M p 



ds 



-1/2 



(6) 



2.1. Flux tube 



Applying the approximations listed in §2 of iBogdan et alJ (jl996l ). a thin flux tube carrying a magnetic 
flux of $>f — 3.88 x 10 17 Mx, with plasma-/3 = 1 everywhere inside the tube is placed in the polytrope. The 
thin flux tube approximation, 

b( s y 



b(s) 



1 8irp(s) 



1 + 



(7) 



where b(s) and R(s) are the magnetic field and the radius of the tube at depth s, is shown to be accurat e 
to better than a percent in the truncated polytrope situated below z = — zq or s = 1 |Bogdan et al.lll996|) . 
Note that the magnetic flux associated with the tube is held constant - different values of (3 therefore result 
in different b(s) and R(s). 
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2.2. Oscillations of the tube: the kink mode 



Horizontal motions of the flux tube created by the impinging modes (whose di splacement potentia l is 
given by ^(x, <)) in the direction of the wave vector are described by £(s, t) (e.g. iBogdan et al.lll996l ), a 
solution to the differential equation 

2(1 + /3) a 3 * 



d 2 



-0 



2.g s 



d 2 



d 



dt 2 (1 + 2(3)(m + 1) ds 2 1 + 2(3 ds 



1 + 2(3 Zo ~dxdt 2 ' 



(8) 



where x = rcosO. Following IBogdan et al 



( 19961) . we define £_L = ~i£(s), where £(s,t) = £(s)e _iwt , and 
£(s) is the purely spatial component of the tube displacement. The function £j_ contains all of the scattering 
information that is needed to understand the interaction of the wave field with the flux tube. 



2.3. Jacket modes and a lower boundary 



Bogdan fe Callv (|l995h showed that scattered waves created as a consequence of magnetic interactions 



are a mixture of modal and evanescent components. The scattering process results in not only a redistribution 
(and loss) of modal energies but also in the production of a continuous spectrum of evanescent 'modes' 
called jacket modes. The existence of evanescent waves in the near field of an arbitrary scatterer (not 
just magnetic) that poss esses subwavelen gth features has been well do cumented in the area of acoustics 
(e.g. Mavnard et al.lll985[ ) and optics (e.g. Bozhevolnvi fc Vohnsen 19961) . Mathematically, these evanescent 
modes appear to compensate for the inability of the set of propagating-mode eigenfunctions to represent 
a spatially intricate scatterer. These non-propagating evanescent waves decay exponentially rapidly with 
distance from the scatterer. 



The formalism of IBogdan fc Callvl (jl995l ) requires an uncountable infinity of jacket modes to complete 
the basis, a consequence of the lower boundary being placed at s = oo. Jacket modes in this instance are 
mathematically described using Whittaker functions, which are relatively diffi c ult an d expensive to compute 
accurately. Moreover, the jacket mode equations listed in IBogdan fe Callvl (|l995l ) contain integrals over 
running indices that pose significant numerical hurdles because of the poor convergence properties of the 
integral. In order to circumvent any calculations involving these continuous jacket modes, we employ a lower 
boundary placed at a depth s = D = 250 (z = 98 Mm) to reduce this uncountably infinite s et of evanescent 
modes to a more tractable discrete countably infinite counterpart in the manner described in lBarnes fc Cally 

(boooh . 



Only modes whose inner turning points are in the vicinity of the lower boundary are affected by its 
presence and therefore, placing it at a depth of 98 Mm means that the / and first few p modes remain 
unaware of its presence (because of the decay of the eigenfunction) . Moreover, the scattering amplitudes 
decrease sharply with increasing radial order. Consequently, the high order p modes which interact with 
the lower boundary are largely unimportant in any case because of their weak contribution to the scattering 
process studied here. With the introduction of the boundary, the problem becomes numerically well defined 
as well. 

In order to obtain the p- and jacket mode eigenfunctions, we apply a zero Lagrangian pressure p ertur- 
bation lower boundary condition for the sake of simplicity (see the appendix in lBarnes fc Callvl 120001 ). We 
also assume that the tube oscillations in Eq. [8]) are oblivious to the lower boundary. This condition is 
necessary because energy loss in the form of propagating Alfven waves along the tube can only occur when 
the lower boundary is transparent (not the case with zero Lagrangian pressure perturbation) . 
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To summarize, we use the formalism and model of Bogdan et al. (199p) but unable to carry out cal- 
culations of the necessary jacket modes described in IBoedan fc Callvl dl995h. we replace this u ncommonly 
difficult continuous set of modes by its more tractable discrete cousin (jBarnes fc Calivll2000h . Next, we 
describe the tube radius boundary condition that allows us to begin the task of estimating the scattering 
coefficients. 



2.4. The scattered wave and tube boundary 



The scattered wave is given by (e.g. IBoedan fc Callvlll995f ): 



E E « 

n— m— — oo 



o i[m0— cut] 



(9) 



where ot v mn arc the p-modc scattering coefficients, f3f nn are the jacket mode coefficients and K m (w) is the 
if-Bessel function of order m and argument w. The resonant wavenumber of the p n mode is denoted by 
kP (including the n = / mode), while wavenumbers corresponding to jacket modes are labeled k 3 . The 
eigenvalues n J n > and K v n are related to k^ and k? respectively, according to equation (|5|) (replace p by J 
in Eq. [5]). fn cquation[3l The un-normalized jacket mode eigenfunction, <&,/(k^;s), is given by 



*j( K i;s)=s-V 2 -» 



M_ 



(10) 



where rj^, a parameter and k„, the jacket mode eigenvalue are determined by the boundary conditions (see 
appendix A for details) and M K)fl (w) is the Whittaker function (jWhittaker fc Watsonlll980l) . The normal to 
the tube boundary at a given depth s is given by 



r - 

zq as 



1 zq as 



1/2- 



(11) 



A consequence of the thin flux tube approximation is that | | < < f ; however we retain this term throughout 
in order to verify that the scattering coefficients are at best weakly dependent on it (see 32]). The boundary 
condition is then obtained by matching radial velocities across the tube boundary, r — R(s), 



lr=R(s) 



x ■ r£±(s)e 



where x is the unit vector along the x-axis. Simplifying equation (|12p . we obtain 



Or 



~v \ 1 dR d , T , ' 

(* inc + * sc ) - ^"^-^(*mc + *sc) 



*£±(s)cos0 



-R{s 



1 dR 

Zq ds 



1/2 



(12) 



(13) 



Using a least squares approach, equation (fill)) is then solved to obtain an estimate for the scattering coeffi- 



cients oP mn . 



3. Solution procedure 



By only retaining the m — 1 coefficients and canceling the e %UJt term in equation (| 13[) . time and angular 
dependencies may be eliminated. Because we are only dealing with a single m, the scattering coefficients 
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introduced in equation ([9]) are rewritten as a^,/3^. We define the following functions 



/(«) = 




1 dR d 



--) + ( g-r ~ ) lMk P Q r)%K;s)} \ r = R 



ds ds 



?.( S ), 



9 P n{s) 



d_ 

dr 



1 dR d 

z% ds ds 



--R(s), 



j id 1 dR d 



dr 



and arrive at the least-squares problem: 



ds ds 



[Ki{k J n r)$j{n J n ;sj\ | r=iI(s) , 



A 




(14) 

(15) 
(16) 

(17) 



In equation (|17[) . si and sm (iVf is the number of grid-points in depth) are the start and end of the finite 
vertical domain, [a], [0\ are column vectors of the scattering coefficients, a p n and (3^ and the Mx (iVi+JVa+l) 
matrix A is given by: 



5o( s i) ■- 9 P Nl (si) 9i(si) 



9n 2 (si) 



A 



,9o(sm) ■■■ g P Nl (sM) 9i(sm) ... 9n 2 {sm), 



(18) 



where N± + 1, N% are the number of p- and jacket modes included in this calculation, respectively. The grid 
spacing in the s space was s et at As = 0.249 or c orrespondingly, Az = 97 km. For physical parameters, we 
utilize the model described in lBogdan et al.l (|19961 ). wherein m = 1.5, $ = 3.88xl0 17 Mx, R(s = 1) = 100 km, 
zq = 392 km, and incident modes of frequency, u> = I-kv^v — 2,3,4,5 mHz, and f3 = 0.1,1,10. The lower 
boundary of the box is set at D = 250 or 98 Mm. The thin flux tube approximation may be invoked to ignore 
the dependence of various terms on derivatives with respect to s in this calculation. Although not shown here, 
we verified the validity of this approximation by demonstrating the invariance of the scattering coefficients 
to the presence of the derivative terms. That the derivative terms contribute very little to the solution is 
an important test of the thin flux tube approximation. Also, we considered a near-field approximation of 
the terms Ji{k^r), Hi(k%r), and Ki{k^r) in equations (fT4|) - (fT6|) . Applying this approximation to J and 
H does not significantly change our solutions. However, the near-field approximation to the Bessel K leads 
to substantially different solutions for the scattering coefficients. 

The coefficients for the first 4 modes were seen to be stable to changes in the depth of the lower boundary, 
as discussed in 32.31 To test the accuracy of the p-mode scattering coefficients, we attempted to match the 
energies of the ingoing with the outgoing and Alfven waves. However, because their amplitudes are so low, 
we find that small changes in the phase of the scattered / mode result in significant changes in the amplitudes 
of the high-order scattered p modes. This implies a lack in the robustness of the result; hence we only show 
the values of the /- and pi-mode scattering coefficients. 

Large numbers of jacket modes (463, 675, 820, 857 for v = 2, 3, 4, 5 mHz respectively) are required to fit 
the right hand side with the accuracy of Figure [T] Least-squares fitting was performed using the backslash 
command in MATLAB. The various complementary and standard Whittaker functions were computed using 
CERNLIB, a freely available suite of mathematical functions. The absolute values of a sample of the 
scattering coefficients, \a^\ and corresponding phases, arg(a^) for plasma-/? = 1 and v — 3 mHz are listed 
in the second and third columns of table [TJ We do not display the jacket mode coefficients here. 
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Our confidence in these values is strengthened by the consistency in the values of the scattering coef- 
ficients over a large number of numerical experiments. We show in the upper panels of Figure [T] that the 
combination of resonant and jacket modes captures the right hand side very well. In the lower panels, the 
contributions of the jacket and resonant modes are separated to illustrate that jacket modes are a non-trivial 
component of the scattering process. Shown in Figure [2] is the agreement with intuitive expectation that the 
higher the radial order of the p mode, the lower the contrib ution. This i s in line with the idea of a diagonally 
dominant scattering matrix (for scattering matrix, see e.g. Braunlfl995l k i.e., a scattering matrix with weak 
off-diagonal terms that decay very rapidly. 

There is a definite dependence of the magnitude of the scattering coefficients on the plasma-/? and the 



frequency of the incident wave (e.g. iBogdan et al.l Il996h . It emerges from our calculations that the wave 



flux tube coupling becomes stronger as the frequency increases (see Figure [2]). Moreover, as the flux tube 
becomes magnetically dominated (/3 = 0.1), the scattering is enhanced, as emphasized by the magnitudes of 
the scattering coefficients. Correspondingly, the flux tube becomes relatively stiff and the tube kink mode is 
more resistant to the buffeting forces of the interacting modes (upper panel of Figure [3]) . The reverse effect 
is observed for (3 = 10, where the tube is hydrodynamically dominated (lowest panel of Figure [3]). Also, the 
scattering coefficients increase non-linearly with the amount of flux (and therefore tube radius) at fixed (3 
and v. Since the thin flux tube approximation requires that the tube radius at the surface be smaller than 
a scale height, we can only study a restricted set of radii (see Figured]). 



4. Conclusions 

We have presented a model of wave scattering of f a magnetic flux concentration. Earlier efforts to model 



magneto-acoustic interactions (e.g. lGizon et al.ll2006f ) have been under less realistic conditions because of the 
related mathematical difficulties. However, it is noted that this model does not capture the complexities of 
the solar case where phenomena such as radiative heat transfer, azimuthal asymmetries in the tube, the lack 
of applicability of the thin tube approximation in near-surface regions, local changes in source properties, 
downflows, etc. probably play an important role in affecting the wave field. Even under relatively idealized 
conditions as in this model, the manipulations required to recover the scattering properties are non-trivial. 
Perhaps this indicates that the formalism is sub optimal when keeping in mind the relatively significant 
effort required to compute in this basis set of Whittaker functions while coping with the restricted generality 



wave interactions with magnetic fields (e.g lCameron et al 



Cameron et al. 


2007: 


HanasoEe 


2007) 



these problems in greater generality. The estimates for the degree of scattering obtained here can be used 
to test the validity of numerical forward calculations. 

It is shown that the contribution to the scattering process by the so-called jacket or evanescent modes 
Table 1. Scattering coefficients for j3 = 1, v = 3 mHz for a flux tube with radius 100 km at s = 1. 



Mode Amplitude Phase 



/ 0.0117 49.6° 

Pi 0.0007 76.7° 
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Fig. 1. — The upper panels show comparisons of real and imaginary parts of the fit (labeled 'matched' in 
the figure) to the corresponding real and imaginary parts of the right hand side (= f(s), Eq. [14]) for /3 = 1 
and v = 3 mHz. The curves are indistinguishable. Contributions to the real and imaginary parts of the right 
hand side (RHS) respectively from the jacket and resonant modes are shown in the lower panels. The lower 
boundary was placed at a depth of 98 Mm; the surface gravity mode with 8 normal and 675 jacket modes 
were used in this matching. The scattered amplitude is normalized with respect to the incident /-mode 
amplitude. 
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Mode scattering amplitudes 



0.025 




Frequency, v [mHz] 



Fig. 2. — The amplitudes of the scattered / (solid line) and pi (dashed line) modes at various values of 
plasma-/? and frequency, v. The p\ scattering coefficients are seen to be substantially smaller than the 
corresponding ones for the / mode. 
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Fig. 3. — Kink modes at v = 5 mHz for different values of j3. The magnetically dominated case with (3 = 0.1 
is seen to be very stiff, as opposed to the (3 = 10 case. The amplitude of the mode is normalized by the 
amplitude of the incoming / mode. 
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Fig. 4. — Scattering coefficients for v = 3 mHz, [3 = 1 as a function of tube radius. The flux, $ oc R 2 , 
where R is the tube radius. A clear trend with increasing tube radius (or flux) is seen. The thin flux tube 
approximation requires the tube radius to be smaller than a scale height (= 260 km); therefore, we cannot 
study the scattering at large values of R. 
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is entirely non trivial, and therefore inevitable in st udies of scattere r s that contain subwavelength variations 
in their spatial structure. It has been realized by iMavnard et al. that these near-field evanescent 



modes may be analyzed to image with subwavelength resolution. Thus, future prospects of studying these 
flux tubes in greater detail are exciting, especially in the context of upcoming satelite missions. 



Lastly, to say something about the nature of the average flux tube in the measurements of iDuvall et al 



(|2006l ) based on this analysis, we first assume that the interactions are well captured by the single scattering 
approximation. This implies that the scattering magnitude will scale linearly with the number of model 
flux tubes while the phase of the scattered wave is unchanged from the value listed in Table [TJ While our 
scattering model over estimates the amplitude of scatter by a factor of 8.8, the phase of the /-mode scattering 
coefficient agrees well with the observations. Several possibilities arise: either the (3 3> 1 for flux tubes in the 
Sun, or the upper boundary conditions of this model are inaccurate, or there are important pieces of physics, 
like non-linearities, radiative heat transfer etc., not accounted for here. More conclusive results await further 
investigations. 

S.M.H. was funded by NASA grants MDI NNG05GH14G and HMI NAS5-02139. Computations were 
performed on the Stanford Solar group machines. 



A. Eigenvalues 

A.l. Jacket mode eigenvalues 

The functional form of jacket modes used in our calculations is given by equation (|10|) . These functions 
are forced to satisfy the boundary conditions, 



<9$ 
~~ds 

at s = 1 and s — D. Defining 



J 



— $j = (Al) 
m 



Nj(K J n ,»,s) = y^LM^^ l^js)- ^M'_ iK] [^js\- ^-M_ lKi ^ l^js), (A2) 





%s) 


iv 2 




K*n J 


K n 



n k„ / m \ k 
and 

r-, / 7 x i/S + M,, (iv 2 \ iv 2 ^, (iv 2 \ v 2 „ r (iv 2 \ ... 

Dj( K J n ,n,s) = -^-Jt M _ i<tll i^-jsj + ^jM'_ i<ttl [^jsj + ^aj , (A3) 

the eigenvalues are determined through the relation 

Nj(Ki, f i,l)Dj( K J n ,^D)^Nj(4, f i 1 D)Dj(K J n ,^l). (A4) 
Subsequently, the constant r\ J n in equation (|10[) is obtained: 

j _ Nj{ K £,n,l) _ Nj{K J n ,n,D) 



A. 2. p-mode eigenvalues 

The functional form of p modes, given by equation ([4|, has to satisfy 

3% | v 2 
ds m 



(A5) 



% = 0, (A6) 
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at s = 1, D. Following the formalism of appendix lA.il we define N p , D p as: 



N p (k,p, h, s) 



1/2 + /i 



11/ \ 1/ . 11/ \ 1/ 11/ 



and 



1/2 + A t 



(A7) 



(A8) 



and determine the eigenvalue K v n and constant CJ^ in equation ((4J through the following relations, respectively: 

N p { K 1,^)D p {kI,^D)^N p {kI,^D)D p {kI,pl,1), (A9) 

CP = a^«,m) = JVpK,M,£>) (A10) 

1) D p («;^,/x,D) 



B. Determining the amplitude and phase of the scatter from observations 



Through analyse s of Michelson Dopple r Imager (MDI) data (jScherrer et al.lll995l ) and the application of 



the forward models of lGizon fe Birchl (|2002( ), lDuvall et al.l (|2006f ) estimated the complex scattering amplitude 
(d in their paper) due to network magnetic magnetic elements. Assume an incident plane wave of the form, 
i.e. an alternate form equation (|3|), 

$ in =e ik ° x - iut %{kl;s), (Bl) 
where x — x ■ r, with the velocity of the wave given by, 



vr = — V$ m = ujk $' 

h dt 



(B2) 



(B3) 



For dipole scattering only, the scattered / mode takes on the form: 

$ sc = -2ial Q cos 9 Hi {k r)%(kl;s). 

Using conventions from Gizon fc Birch ( 20021) and the definition of dipole scattering of Duvall et al. ( 2006 ) , 
we have, 

C Sv z = -6£v in = -V h • [dt#] . 
Solving equation (|B4[) for a point scatterer, 

Sv z (r, w) = -2vr / ds G z (r - s, w) V s • [dS(s) <(s, w)] , 

which upon integrating by parts produces 

5v z {r,uj) = 2nd / ds 5(s)< 1 (s, w) • V s G^(r - s, w), 



(B4) 



(B5) 



or, 



5w z (r, w) = -2nd <(0, w) • V r G z (r, w). 
In the far field, the Green's function reduces to 



(B6) 
(B7) 



(B8) 
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The derivative in equation (|B7|) becomes 



V r G z (r)*-^Hl{k r)r, (B9) 
Airg 



and equation (|B7|) along with the expression for the wave velocity from equation (|B2[) simplifies to 

6v,(r,u>) = y p (k p ; S ) l -^H 1 1 (k a r)cose, (BIO) 

and therefore, 

kd 

$ sc = 5$ = 6v z /(-iku) = cos6Hl(kor)<f> p (kP; s). (Bll) 

From equation JB3|, we obtain the following scattering coefficient: 



a P w = - l A d (B12) 



4.9 



For d with a phase of about 130 deg ([Duvall et al.ll2006f ). a\ has a phase of about 40 deg 
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